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1 Introduction 



The quantum chromodynamics (QCD) is a theory of the strong interaction. So far 
all experimental data are consistent with the predictions of QCD. Especially the 
high energy behavior of QCD is believed to be described by the perturbation theory 
thanks to the asymptotically free nature of QCD. Among many interesting and 
important problems, the spin structure of nucleon has been one of the most exciting 
subjects in recent years. The data on the polarized deep inelastic scattering by 
the EMC collaboration result in excitement of not only particle physicists but 
also nuclear physicists since the data seem to indicate that the nucleon's spin is not 
carried by quarks (partons). This EMC experiment has incited many (particle and 
nuclear) physicists to challenge the so-called "spin crisis" problem in QCD. After 
a flood of theoretical papers as well as new experiments our understanding 

on this problem is now much more improved: the interpretation of the QCD results 
in terms of the parton model is never obvious and simpleminded one may fail: the 
axial anomaly plays an important role: etc. 

The aim of this talk is to provide a pedagogical introduction to the perturbative 
QCD to study the nucleon structure through the deep inelastic process for non- 
experts of QCD who are interested in the deep structure of the nucleon. Those who 
are familiar with QCD and interested in recent progress in this field are referred to 
recent nice article and reviews [0. 

In Sect. 2 we review the kinematics of the deep inelastic lepton nucleon scattering 
process. In Sec. 3 the basic approach of the perturbative QCD based on the operator 
product expansion and the renormalization group equation will be explained. The 
relation between the QCD results and the parton model is discussed. It will be 
stressed that the parton density is a "conception" which depends on the renormal- 
ization scheme. In Sec. 4 we consider the polarized structure functions. Concluding 
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remarks including some subtleties and/or controversial aspects which deserve more 
investigations to understand recent experimental data are given in Sec. 5. 

2 The structure functions 

The cross section for the deep inelastic lepton {l{k)) nucleon {N (p)) scattering 
l{k) + N (p) I (k') + X (Fig.l) is given in terms of the leptonic and the hadronic 

X 



1 N 

Figure 1: Kinematic variables in inelastic lepton- nucleon scattering. 

tensors according to the standard procedure in the field theory. 

da _ 1 f 
^d?k' k-p\47rQ^ 

where we consider only the QED interaction between the lepton and nucleon and 
keep only the lowest order in agED- ? is the momentum transfer from the lepton to 
the nucleon and = —Q^ = {k — k'Y . The leptonic {L^i,) and the hadronic (W^jy) 
tensors are defined as follows; 

V = \ Eik, s\UO)\k', s') {k', s'\UO)\k, s) , 

s' 

W,. = ^Y.{P^S\J,mX){X\JM\P,S){27r)'6\px-p-q) 

ATT ^ 

= ^Jd'xe^^-^{p,S\[J,{x),MO)]\p,S) . (1) 
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with the lepton's (hadron's) electromagnetic current {Ju)- s{S) is the spin of the 
lepton (nucleon). 

In general, the spin 4-vector of the fermion with mass m and momentum k is 
defined as s ^ = 1 , = — 1 , s ■ A; = 0. So, for the longitudinally polarized (helicity 
± ) states, we get, 

s^ = ±— (A;, 0, 0, A;°) , A; = I^P . 
m 

We can use the following approximation for leptons ms'* ~ since miepton — 0. 
Using this approximation, we get for the leptonic tensor, 

± I I (f X a 

^IJ,v ~ ^fJ^v "I" ^vk^ + -^g^iv -F 'i'^ixuXaQ k'^ ■ 

On the other hand, the hadronic tensor contains all the information of the strong 
interaction (QCD). Taking into account the various symmetries, namely the Lorentz 
invariance, current conservation of the QED current, T and P invariance, we can 
write down the general form for Vl^^,^, 



with 



( .j,r , ( p-q \ [ p-q \W2 

= s,.x.q^[s''MG, + (p-qS''-q-Spn^y 



where M is the mass of the nucleon. W^^ i^j^u) is the symmetric (antisymmetric) 
part in /lu and relevant to the unpolarized (polarized) process as shown below. Usu- 
ally we define the following dimensionless structure functions ( scaling functions ) : 

u Mu z/2 

with Mv = p ■ q. These structure functions depend on the Q'^ and u or Q'^ and x 
(Bjorken variable) x = ^ — ^j^- 
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The explicit formula for the cross sections for this process in the Laboratory 
frame corresponding to the configuration in Fig. 2 is easily calculated to be 0, 




Figure 2: Momentum and spin configuration in Lab. frame. 



da^'^ da _^ da^ 



with 



dE'dn dE'dn dE'dQ ' 



- ' 2iyism - + 1^2 cos^ 



dE'dn Q^M \ 2 2 

da"^ a^E' 



cosa {{E + E' cosO) MGi - Q^G2} 



dE'dn MQ'^E 

+ sin a cos E' sin e{MGi + 2EG2] 

The superscript ± refers to the lepton's helicity and E' = k'^ and E = k^. From 
these formulae, we can derive the expressions , for example, for the longitudinal 
asymmetry which is the difference between the cross section for the nucleon's spin 
being parallel to the lepton's (H) and the nucleon's spin being anti-parallel to the 
lepton's (Ti) : 

da^^ da^^ 2a^E' 



{mGi{E + E' cosO) - Q^G2} ■ 



dE'dn dE'dn MQ'^E 

It is traditional and sometimes convenient to express the structure functions in 
terms of the virtual photoabsorption cross sections (in Lab. frame). The definition 



of the virtual photoabsorption cross section is given by, 

7r6^ 

with the photon's polarization vector e^. We have adopted the Hand-Berkelman's 
convention. Taking the direction of photon's momentum q to be z-axis ( note that 

— * — * 

q k ), we have the following polarization vector for photons in the nucleon's rest 
frame: 

q^ = (u, 0, 0, ^Jly' + Q') , 

I , . 1 



en = -^(0, 1, 0) , £5 = ^(V^' + 0' 0' ^) • 

For the unpolarized structure functions Wi and we get the relations: 

,2 \ 



Tre^ 



(75 = 



2M{u-Q^/2M) 



Gt — W^. 

^ 2M(i/-gV2M) ' 

where T — R and/or L. The unpolarized lepton-nucleon scattering cross section is 
given in terms of as and cry, 

where 

Here we note that e means the ratio of T- and 5'-photon present in the virtual 
photon . 

To get the expressions for the polarized structure functions Gi and G2, let us 

consider the three types of the photoabsorption processes: (T1/2 (cl with S'^ = 
(0, 0, 0, ±1)); (e« with 5''^ = (0, 0, 0, ±1)); (Tts (the interference between T- and 
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S- photon with Sy = ±1). It is easy to obtain, 

2 

Note that, ar — |(ci/2 + C3/2) • By defining the asymmetries Ai and A2: 
_ (Ti/2 - (T3/2 _ MuGi - Q'^G2 ^ _crTS _ MGi + PG2 



ci/2 + (J3/2 W^i cft Wi 

we can write the longitudinal asymmetry A as; 



where 

1 + sR ' ^ E - E'e ' gt 

D is called as the depolarization factor and its physical meaning is obvious in 
Fig. 3. It is easily verified that. 



> Z 



Figure 3: Momentum configuration in inelastic lepton-nucleon scattering. 
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with 

1 - {E'/E)e 



Di = cos ij) 



= probability to have T— photon — 

gt + sag 

The structure functions gi and 5^2 can be expressed in terms of Ai and A2. 

«' = 2j(rb) + • * = 2j(rTfl) (^•''-^ • 

where the following relation has been used, 

;(if«) 4) ■ 

Here it is to be noted that we have the inequalities from the unitarity arguments ; 
1^1 1 < 1 , 1 < \/R. Furthermore if we consider the scaling region (Bjorken limit), 
^ = 1^ « 1 , is given by. 
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2x{l + R) 

This is the basic formula on which the experimental determination of gi is based. 



3 The perturbative QCD 

In this section, we review the fundamental aspects of QCD to analyze the structure 
functions introduced above. At first, we discuss the general strategy of QCD based 
on the operator product expansion (OPE) and the renormalization group equation 
(RGE). Next, we consider the relation between the QCD results and the parton 

model interpretations of the process we are considering. 
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3.1 Formal approach in the perturbative QCD 

The kinematical region which we are interested in is the Bjorken hmit, namely 
, z/ ^ oo with X = 2^ fixed. In this hmit, we can easily recognize that the 
hadronic tensor Eq.(|l|) is governed by the behavior of the current products near 
the light-cone. So the light-cone expansion, which is a variant of the OPE, of two 
currents might be applied to this precess. However the OPE can not be used directly 
here by the following reason. The OPE makes sense in the short distance limit and 
this limit corresponds to the region where all component of become infinity. 
Therefore, >2p- q. On the other hand, the physical region for the deep-inelastic 
process is < 2p ■ q. Fortunately we can overcome this dilemma by using the 
dispersion relation |Q] which relates the short distance limit to the Bjorken limit for 
the deep-inelastic process. 

Consider the time-ordered product of two currents which corresponds to the 
forward Compton scattering of the virtual photon with "mass " , (the Lorentz 
and the spin structure being neglected), 

T = iJ dxe''^-''{p\TJ{x)J{0)\p), 

the physical region of which is 2p ■ q/\q'^\ = 2Mv/Q'^ < 1. For this process, we can 
use the OPE. In general, the product of two operators (currents) can be expanded 
as follows; 

Tj{x)j{o) ~ E^r(^' - ^^)^"^^"^ ■ ■ ■^""o;^...,„(o), (2) 

i,n 

where C" is a c-number function called Wilson's coefficient function and O* ,, „ 
are local composite operators labeled by the index i. The dimensional argument 
tells us that C"(a;^) behaves like. 
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in the concerned limit of ^ with x'^ small but 7^ 0, where d is the dimension 
of the corresponding operator. So the operators with the lower twist {jn): 

Tjv = dim. — spin = (Iq — n 

dominate. The Fourier transform of Eq.(^ assumes the following form with an 
appropriate normalization to define Cf (Q^), 

By defining the matrix element of O* , 

(plC'Mi.-Mn(O)b) = 2^>Mi ■■■P^^n- trace terms, 

(note that the operators should have definite twists, so must be traceless) the forward 
Compton scattering amplitude T is written as, 

i,n ^ ^ 

To get information for the structure functions in the Bjorken limit, we rely on 
the analytic structure of T in the complex tu-plane. There are cuts going out to 
infinity from uj = ±1 and the discontinuity of T is related to the structure function 
W: W = ^ ImT. The next step is just to use the Cauchy's theorem. The contour 
integral of T around the origin in u;-plane picks up its n-th coefficient: 

Deforming the contour to pick up the discontinuity of T, the left hand side becomes, 

LHS = - / dujuj-''-HmT{uj,Q^) = 2 / dxx''-^W{x,Q^), 
n Ji Jo 

where we have used the crossing symmetry, 

T{q,p) = T{-q,p) , W{q,p) = - W{-q,p). 

9 



Finally we get the so-called moment sum rule, 

Jo Y 

(In the above derivation, it is assumed that the integral is convergent at infinity. 
Since this region corresponds to the Regge limit, the problem of convergence is 
controlled by the Regge behavior of the corresponding amplitude.) 

Now the Wilson's coefficient functions C^iCp) depend on only the large momen- 
tum Q^. So the RGE is effectively apphed to Cf (Q^) and they are calculated in the 
perturbation theory due to the asymptotic freedom of QCD. On the other hand, the 
hadronic matrix element of the composite operator can not be estimated pertur- 
batively (long distance physics) and treated as input parameters in the perturbative 
QCD. The RGE for the coefficient functions is easily derived as follows. Consider 
the Green's function of currents and fundamental fields 0; 

G)j = (O|TJ(a:)J(O)0(a;i)---0(a;,)|O) 

= ^C"(x)(O|rO„(O)0(xi).--0(x,)|O) = ^C"G^, 

n n 

where we write the OPE symbolically as J(x)J(0) = En C'"(x)O„(0) . The RGE 
for LHS reads, 

where 7j , 7 the anomalous dimensions of J , . P is the well-known operator in 
the obvious notation, 

with m the quark mass. The RGE for RHS becomes. 
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with 7q the anomalous dimension of 0„ . 

Therefore RGE for the coefficient functions is given by, 

[P + 27j-7a C" = 0. 

It is easy to solve the RGE and the solution for the Fourier transformed coefficient 
function becomes, 

C"(Q^) = C'^(^^,gii,),m{i,)^^C"il,m,fn{t)) 

xexp fdt'[2^j{g{t'))-YMt'))], 
Jo 

where, 

i = (l/2)ln(gV/.^), 

dg dm /-xx - 

with (?(0) = (^(/i) , m(0) = m(/x). /i is the renormalization point. It is to be noted 
that the above solution takes, in general, a matrix form because the operators with 
the same quantum number and twist will mix under the renormalization. 

Here let us consider how to calculate the anomalous dimension 7 and the "coeffi- 
cient function" , C"^(l, ^(t), m(t)). The anomalous dimension 7 can be obtained from 
the renormalization constant for the composite operator On- To obtain C"( 1, ^, m) , 
we use the fact that the OPE is an operator relation. So, consider a Green's function 
which can be explicitly calculated, for example, 

r(p,g) = iy"(iW'^-^(o|T0(-p)j(x)j(o)0t(p)|o) 
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The LHS, the current correlation function, will be calculated to be, 
r(p,g) = 2^t„u;" , tn = tl + gX + ---. 

n 

The RHS becomes, 

n 

where a„ is obtained from 

(O|T0(-p)O„0"^(p)|O) = 2a„p^^ ■ ■ -p^^ - trace term. 

Therefore, 

tn = Ctv.C"' . (3) 

Expand C^iQ"^) in powers of 



9^ 



2 



9=0 



where 7q = 79(7^ + 0{g^) and, for simplicity, we neglect the mass m and take into 
account the fact that the anomalous dimension of the conserved current J vanishes. 
By comparing terms of the same power of g in Eq.(|^), we will get C"( l,g). 

Now the realistic case of QCD. The OPE of the electromagnetic currents looks 
like. 



X 
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The explicit forms of the composite operators with the lowest twist which contribute 



to the unpolarized structure functions |T0[ read : 
O 



o 



G 



where A* are the SU{f) generators, D is the gauge covariant derivative, and S 
means the symmetrization on the Lorentz indices. To the spin dependent struc- 
ture functions, the following twist-2 and the twist-3 {R2) [|12| operators 



contribute [|T3|. 





■■fJ.n-1 




4ay, 




td1,F 








(4) 













.2,F 

'0"Ml---Mri-l 



n 



[(n - l)V'757<x^{/.i ■ ■ ■ 

n-l 

5^ ^757/.i^{a^Mi ■ ■ ■ ^W-i^rt+i • ■ • ^Mn-l}V^ 



(5) 



2,m 
2,fc 



where 



V'fc = i'"#75^/.i ■ ■ ■ G^^., ■ ■ ■ /^^„_27Mn-iV' , 
f/fc = I'^-^g^D^, ■ ■ ■ G„^, ■ ■ ■ D^„_,7Mn-i • 

We have shown only the flavor non-singlet operators for R2. 
Defining the matrix element of composite operators by. 
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-IJ-n- 


Jp,s) = 


-2MaiS{„Pf,^---p^^_,} 


(P,S\RZ- 


■■Hn- 


Jp,s) = 


71/ — 1 


{p,s\Rl';r,. 




Jp,s) = 


-2 Men{S^Pi_,, - S^^pa)p^^ ■ ■ ■Pi_,^_, , 


{P,S\Rli. 


■■/J-n- 


Jp,s) = 


-2 Mf^{S^p^, - S^,p^)p^, ■ ■ ■p^^_, , 



The moment sum rules are written as, 

3 3 

j dxx^-'g,{x,Q') = \ Y.<EI^{Q')^ 
J dxx''-^g2{x,Q^) = 



3 

n — 1 



1 



+ - 



dtY{0)) 



k 

The each coefficient function ( E"^ or ) takes the following form neglecting the 
operator mixing, 

E-{Q')=E-{lrg{t))ew 
If we insert the perturbative results, 

rig) 

Pig) = -Pog'-Pig' + 

we can get the final answer at the one-loop level, 
where 



7o / + 7ig' + 



(6) 



2/?o Wl ) j 

-7? /2/3i +7^/2/30 



g 



7o"//3o 



(7) 



In this way, the dependence of the structure functions can be predicted in the 
perturbative QCD based on the OPE. 
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3.2 Parton pictures in the perturbative QCD 

It is now a widespread belief that the resuhs of the perturbative QCD (at least at the 
lowest twist level) can be understood in terms of the "parton "language. In fact, the 



idea of the QCD improved parton model ||Tj] has been justified for various processes 



and has produced a great deal of progresses in many hard reactions, especially those 
to which we can not apply the OPE directly. As far as the deep inelastic process, 
which we are interested in, is concerned, we do not need the parton model at all 
since we can make a definite prediction based only on the OPE and RGE. To relate 
the QCD results based on the OPE with the parton picture, we must define the 
parton distribution function in a appropriate way . 



Let us define the parton distribution function from the previous results based 
on the formal approach of QCD. Consider, for example, the moment of F2 in the 
singlet channel. 



The QCD says. 



where 



Texpj- dt'^-{-g{t'))] 



and An{Q'^) is the nucleon matrix element of the composite operator renormalized at 
Ql- The Ql dependence is (should be) cancelled between An{Ql) and the coefficient 
function. Since Qo is arbitrary, let us put Ql = Q^- 

M^iQ')=Af^{Q')CUl,-9iQ'))+A^{Q')C^il,giQ')) . 

Now we can define the "Q^ dependent parton distribution function" of the quark and 
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gluon; q{x,Q'^) , gi^x^Q"^) as, 

A^{Q^) = I dxx''-\{x,Q^) , A^(g') = j dxx''-'gix,Q^) . 

Expanding C" perturbatively , 

C^lrg) = 1 + clj^f + ■ ■ ■ , = claf + • • • , 

and making the inverse Mellin transformations, we get, 

-F,{x,Q') = qix,Q')+f ['^KJ^]qiy,Q') 
X Jx y \y ) 

+f f-Koi-]9{y,Q') + 0{f\ 
y \y/ 

where 

J dxx''-^KF{x) = dip , j dxx^'-^Kcix) = c\q. 

We can interpret the above expressions as: the first term describes the interaction of 
the charged parton (quark) at the lowest order and other terms indicate the radiative 
corrections for quarks and gluons at the higher orders. So, the QCD results can be 
interpreted in terms of the parton model with the dependent parton distribution 
functions. Furthermore we must take into account the presence of the radiative 
corrections. 

It is here to be noticed that the above definition of the parton distribution 
function (consequently, the corresponding coefficient functions too) has a degree of 
arbitrariness. Since the parton distribution functions are defined as the nucleon ma- 
trix element of the composite operator renormalized at , they depend on how to 
renormalize the composite operators. To understand the source of the arbitrariness 
in the definition of parton distribution functions, let us remember the procedure of 
obtaining the coefficient functions. In general, the quantities in Eq.(^) will be given 



16 



as, 

an = l + ^7'|&„+Q7o"-7F)ln^| + 0(^?^), 

= l + /(ci-^7o"ln^)+0(^7^), 

where we assume = a^ = Cq = 1. is the anomalous dimension of the external 
field with which one estimates the both sides of the OPE. So, from Eq.(|^) we get; 
c" = dn — hn- Now hn depends on the renormalization scheme adopted. Therefore 
also depends on the scheme. The scheme ambiguity comes from how to renormalize 
the composite operator. Suppose that the bare expression for a„ to be. 

To renormalize a„, we multiply ^2^0^; 

where j discriminate various renormalization schemes. Therefore, 

h = f - 

Jn 1 

and this means, 

= dn — {fn — Zl). 

The ambiguity in Z2 cancels between t„ and a„, however 2;^ is quite arbitrary. So, 
the composite operators and also the coefficient functions turn out to depend the 
renormalization scheme. 
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Physical quantities like moments should not depend on the renormalization 
scheme adopted. How does the above ambiguity cancel in the final expression ? 
The difference between the schemes for the composite operator; 



'0„ ' 



affects the anomalous dimension of composite operator. 

d 



7 



n(j) 



= /i— In 

ail 



-2(4-^^)/?o/ + 0(/) + 7 



n{k) 



Then, 



7r^^'^ = -2(4-4)/3o + 7r 



n{k) 



Therefore the ambiguity cancels in the final answer through the combination, 

in C"*((5^). But each term can depend on the renormalization scheme. We can say 
schematically. 



LHS does not depend on the scheme. However the arbitrariness remains in the 
definition of the parton distribution functions. The point is that the use of the parton 
language is helpful to interpret the QCD results intuitively and/or economically. So, 
the definition of parton distribution functions can depend on one's convention. It 
should be also mentioned that even if we start from the parton model without using 
the OPE, the same ambiguities appear at the stage of the factorization. 
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4 The spin structure functions 



In this section, we will give a little bit more detailed analyses for the spin structure 
functions gi and g2. We will not mention to the numerical studies of the recent 
experimental data [1^, l7|. Let us remember the moment sum rules for gi and g2 



explained in sec. 3 [|T^, ^ 

J dxx--'g^{x,Q') = ^-J2<El^iQ'). (8) 
j dxx^-'g2{x,Q') = -"^[al^El.iQ') - dnEl^iQ')] (9) 



1 

+ 2 



k 



The first moments n = 1 of the above equations are related to the interesting sum 



rules: Bjorken sum rule I^TKthe flavor nonsinglet part of gi), Ellis- Jaffe sum rule [22 
(the flavor singlet part of gi) and Burkhardt-Cottingham sum rule ((72)- 

At flrst let us consider g2. For general n, the QCD analysis is not so straight- 
forward for g2 as for gi. The fact that the twist-3 operators also contribute to g2 in 
the leading order of l/Q"^ produces new aspects which do not appear in the analy- 
ses of other structure functions. The appearance of composite operators which are 
proportional to the equation of motion makes the operator mixing problems rather 
complicated. This is a general feature of the higher twist operators. This problem 
has been discussed in Ref. ||20[| and references therein. However, as far as the flrst 
moment {n = 1) of g2 is concerned, the above complexity is irrelevant since there is 
no operators corresponding to n = 1 (see Eq.(||)). So, eq.(P) for n = 1 predicts the 
Burkhardt-Cottingham sum rule: 

J dxg2ix,Q'^) = 0. 

Recently the validity of this sum rule was questioned and explicit calculations of 



g2 (current correlation functions) at the one-loop level have been performed |26 
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The results confirm that the Burkhardt-Cottingham sum rule does not receive any 
radiative corrections in the (at least) perturbative QCD. 

Next we turn to gi Eq.®. Three operators in Eq.(^) contribute to the moment 
of Qi. For convenience of explanation, let us rewrite the operators and their proton's 
matrix elements in the following way: 



R. 



1,0 



R. 



R. 



1,G 



R^ 



1,G 

""Ml-'-Mn-l 



(10) 

(11) 



where i is the the flavor index (we consider the case of 3 flavor i = u, d, s) 



(p, S\Rl;^\p, S) ^ -2 M/\iS{^p,, ■ ■ -p.^^.y. 

Here j runs over i, and G. In this notation, the moment sum rule for gi Eq.( 

reads, 

r(Q2) = ldxx''-'g,{x,Q') 



(12) 



after taking into account the quark charge factor. A° = A'^ + + and (e^) = 
J2i ^11 f — 2/9 for the number of flavors being / = 3. E'^g = E^- is the coefficient 
function for the flavor non-singlet operator Eq. ([TI1|) and does not depend on i. On 
the other hand Eg = E^ q and Eq = E"q correspond to the flavor singlet channel 
Eq ■ ([TI|) and get mixed under the renormalization. 

The QCD calculation of the coefficient functions at one loop level has been 
done many years ago for both the non-singlet [1^, |2^ and singlet [1^ parts. In the 
following, we restrict our discussions only to the first (n = 1) moment. For a nucleon 
target (p: proton, n: neutron), the first moment becomes from Eq.(|T2|), 



12 



-)(A^-A^) + -(A!; + A?-2A^,) 



(13) 
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Since the gluon operator (corresponding coefficient function) with n = 1 does not 
exist (see Eq.(^), only the fermion bihnear operators contribute which turn out to be 
the axial vector currents. Note that the anomalous dimension of the non-singlet axial 
current vanishes because of the current conservation in (massless) QCD. Therefore 
the nucleon's matrix elements of the non-singlet axial currents are scale independent 
and given by, 

A^-A^ = ^ = 1.26, 
+ - 2A\ ^ 0.58 , 

using the flavor SU(3) symmetry. The corresponding coefficient function has been 



calculated in the perturbative QCD up to the three loop in MS scheme |30 



IT V TT y V TT / 

where as = g"^ {Q"^) / An . Since the Bjorken sum rule receives the contribution only 
from the non-singlet channel, we can make a definite QCD prediction for this sum 
rule in the sense that the matrix element is known. 



ri _ ri ^ 1 Ga 



l-^-3.58r^V-20.2 (^]\ 



'V 

The flavor singlet axial current is not conserved due to the Adler-Bell-Jackiw 
anomaly [^. This fact makes the matrix element scale dependent [|l^, |2^. The 
perturbative result for E^Q'^) is given as (see Eqs.(|,^), 

EUir9{t)) = i- — + ---. 

TT 

Due to the anomaly, the anomalous dimension 7 starts at the two loop |T9| , 

7o = ' 7I = (Y^24C2(i?)T(i?). 
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(For the higher order corrections, see Refs.|pi||.) Using these resuhs, we write the 



singlet parts as follows. At first, let us define AS(Q^) by (see section 3.2) 

El{Q')A',{^^)^Elilr9{mJ:{Q'). 

Then, 



AE(Q2) = exp 



'dt'^\g{t')) 







AS(^- 



2^ 



The final expression becomes. 



In the remainder of this section, we will mention to the so called "spin crisis" or 
"spin deficit" [jl^ problem. At first we want to stress that this is never the problem of 
QCD. From the viewpoint of the formal approach, we have already finished all tasks 
to predict the moment of the structure functions. What we should do next is just to 
compare the above results with the experimental data. On the other hand, it is also 
true that the interpretation of QCD in terms of the parton language is very helpful 
and convenient to understand the hard processes. As explained in section 3.2, the 
definition of parton densities depends on the scheme; then is ambiguous. In most 
cases like the unpolarized structure functions, however, this ambiguity produces 
only small differences and a naive parton model interpretation holds rather well (of 
course, with the parton densities). If it is the case also for gi, we will have. 



AS = y dxAqi 
, Jo 



X 



where 

Aqi{x) = qi+{x) + qi+{x) - qi-{x) - qi-{x), 
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with qi± {qi±) being the density of flavor i quark (antiquark) parton with hehcities 
±1/2. So if the nucleon's spin is carried by charged partons, we will expect AE ~ 1. 

The first EMC resuh [Q was AS ~ 0.1 (A^^ 0.2)!! (for the recent status of data, 

see Ref.p|.) 



A resolution to this problem was proposed in Ref. [^. In the approach of the 
QCD improved parton model, they have derived the relation, 

AE(g2) = ^ dx (Ag,(x, Q2) _ /gA<?(x, Q^)) , 

with 

g± is the gluon density with helicities ±1. This result says that AS does not measure 
the spin fraction carried by quarks. Furthermore, it was also shown that the second 
term is independent of the scale in the leading logarithmic approximation. (There 
may be some nonperturbative corrections to the above relation. See Ref. for the 
details.) So if the second term assumes some finite value, a smallness of AS is not 
necessarily in contradiction with our naive expectation. However the story is not so 
simple because we can reach a different conclusion even if we work in the the QCD 
improved parton model depending on the regularization of soft singularities which 
turn out to be the scheme dependence pS] . 

Now how does the above result reconcile with the approach based on the OPE ? 
We do not have a gauge invariant operator composed from the gluon field which may 
correspond to Ag. The answer is the following: we just introduce the "Chern-Simons 
current" a.s n = 1 gluonic operator. 

In this case, we can make a finite renormalization of the composite operators between 
the axial current and k^j. And we can get under a particular scheme that the 
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coefficient function Eq which is associated with takes the value, 

TT TT Z 



Therefore, 



with, 

(P, 5| [P}'°1 ^ |P, 5) = -2MA^5', , (P, 5| [A;,]^ |P, S) = -2MA^5„ 

where the subscript R denotes the renormahzed operator. If we identify A^ (A*^) 
to be the moment of Agfi(x) (Ag'(x)), there is no contradiction. In the OPE of 
the electromagnetic currents, only the gauge (BRS) invariant operators can appear. 
The Chern-Simons current, however, is gauge dependent. This gap will be ffiled, if 
one notes the following trick. In the OPE, only the operator R]^^ in fact appears. 
But we can always write as, 

4*.) 

The first (second) term will be identified as quark (gluon) densities. In fact, in the 
above renormalization scheme. 



-.2 



at the lowest order. 

Now ambiguities can come in the above game. Even if we introduce k„, we can 
get different answer depending on how one renormalizes the composite operators. 
Namely we can put an arbitrary number in the front of the Chcrn-Simons current 
by choosing a different scheme. This fact exactly corresponds to the ambiguities 
which appear also in the QCD improved parton model approach mentioned before. 
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In this way, the definition of the parton densities inevitably depends on the scheme. 
So it is a waste of time to make a discussion without specifying the scheme one 
adopted. The point is : Which scheme will be convenient to parametrize the parton 
densities and easy to understand the physics ? And once one defines the parton 
distribution functions, whether one can explain other processes ? (The definition in 



Ref . ||3^ seems to be reasonable in the sense that the quark helicity Ag is conserved 



to all orders since the corresponding operator is conserved current by construction.) 

5 Concluding Remarks 

We have explained, in this review, only the first stage of the approach in the pertur- 
bative QCD to analyze the structure of nucleons. New experiments at CERN [0, 0] 
and SLAC [^, ^ and detailed theoretical investigations have enormously improved 
our understanding on the (spin) structure of nucleon. Now all data seem to be 
consistent with the expectations from the perturbative QCD. It is impressive that 
we can make a precision test of QCD to the accuracy of about 10% Q for the 
spin structure of nucleons. Furthermore, the problem of "spin crisis" became less 
exciting. We can construct a consistent picture in terms of the parton language at 
least qualitatively for the spin structure function [^, [35[] . 



On the other hand, we would like to emphasize that there still remain many 
subtleties and controversial aspects when one tries to make a much more precisely 
quantitative prediction of QCD. We will list up below several points which will need 
more investigations. To compare the experimental data with the QCD prediction, 
we must take into account many theoretical corrections and sometimes make several 
assumptions. Among those problems are : the small x (and large x ) behavior of 
the structure functions theoretical treatments of the "wee" partons including 
the role of the s quark and the iso-spin SU(2) violation the depend ence 
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of the structure functions both in theories and experiments [Q: etc. Since the 
experimental data on gi are taken at rather small values of Q^, the power corrections 
in to the sum rule may give some modifications. The power corrections come from 
the higher twist terms as well as the target mass effects |3^. For the higher twist 
effects, many authors have tried to reveal their theoretical behavior and numerical 
importance in the analyses of data [0. However it has been pointed out recently 
that there is not only phenomenological but also theoretical problem in including the 
higher twist contributions [0]. Due to the renormalon singularities, there appears 
some degree of arbitrariness in the twist expansion. These problems require much 
more elaborate theoretical investigations. Forthcoming precision measurements will 
also provide us with much information and clues to answer the above questions. 

As explained in the text, parton densities depend on the scheme. Therefore to 
confirm the parton model interpretation of the spin structure function, we need in- 
dependent experiments which could measure the parton (quark and gluon) densities. 
If we consider a quite different process (e.g. Drell-Yan process), we encounter other 



spin-dependent structure functions. The experimental and theoretical studies 



of these structure functions must help us to fully understand the spin structure of 
nucleons. Nuclear effects in measuring the neutron structure functions and nuclear 
structure functions are also interesting and important. We need more investiga- 
tions |4^, on these subjects. 

Finally, we hope that various kinds of new experiments and more detailed the- 
oretical investigations will be able to clarify the not only the perturbative but also 
the nonperturbative aspects of QCD related to the Spin Physics. 
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